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ABSTRACT 

We have produced a cleaned map of the Wilkinson Microwave Anisotropy Probe (WMAP) three-year 
data using an improved foreground subtraction technique. We perform an internal linear combination 
(ILC) to subtract the Galactic foreground emission from the temperature fluctuations observed by 
the WMAP. We divide the whole sky into hundreds of pixel groups with similar foreground spectral 
indices over a range of WMAP frequencies, apply the ILC for each group, and obtain a CMB map 
with foreground emission effectively reduced. With the resulting foreground-reduced ILC map (Figure 
|3p, available on-line), we have investigated the known anomalies in CMB maps at large scales, namely 
the low quadrupole (I — 2) power, and the strong alignment and planarity of the quadrupole and the 
octopole (7 = 3). Our estimates are consistent with the previous measurements. The quadrupole and 
the octopole powers measured from our ILC map are 

respectively. The 68% confidence limits are estimated from the ILC simulations and include the cosmic 
variance. The measured quadrupole power is lower than the value expected in the concordance ACDM 
model (1250 /iK 2 ), in which the probability of finding a quadrupole power lower than the measured 
value is 5.7%. We have confirmed that the quadrupole and the octopole are strongly aligned with 

angle 6*23 = 11?8 + ^, and are planar with high planarity parameters t = 0.981q'q2 f° r I = 2 and 

t = 0.91^0 03 for I = 3. The observed angular separation #23 is marginally statistically significant 
because the probability of finding the angular separation as low as the observed value is 4.3%. However, 
the observed planarity is not statistically significant. The probability of observing such a planarity 
as high as the measured t values is over 18%. The ILC simulations show that the residual foreground 
emission in the ILC map does not affect the estimated values significantly. The large scale modes 
(I = 2-8) of SILC400 shows anti-correlation with the Galactic foreground emission on the southern 
hemisphere. It is not clear whether such anti-correlation occurs due to the residual Galactic emission 
or by chance. 

Subject headings: cosmic microwave background — cosmology: observations — methods: numerical 



1. INTRODUCTION 

The Wilkinson Mic rowave Anisotropy Probe (WMAP; 
iBennett, et al.l l2003al) has measured the cosmic mi- 
crowave background (CMB) temperature anisotropy and 
polarization with high resolution and sensitivity and 
opened a new window to precision cosmology. The 
WMAP one-year data imply that the observed CMB 
fluctuations are consistent with predictions of the con- 
cordance ACDM model with scale-invariant adiabatic 
fluctuations ge nerate d during the inflat ionary epoch 
llHinshaw et all l2003t iKogut et all 120031 ISpergel et all 



120031 iPage et alJ 12003 iPeiris et alJ I2003[K The re- 
cent release of the WMAP three-year data has con- 
firmed the primary results from the one-year data, giv- 
ing more accurate determination of cosmological parame- 
ters llHinshaw et al.l2006l ISpergel et all2006HPage et al.l 
2006). 

However, some peculiar aspects of the CMB maps have 
been noticed, leading to many controversies. First, there 
have been many reports of detection of non-Gaussian 
signatures (iChiang et alJl2003t lParkll2004l lEri ksen et al " 



l2004bl(ilColes et al.H2004HCopi et al]l2004HVielva etal 



12001 ICruz et a,U 120051 120061 iToieiro et all 120(161). as 

oppo s ed to the WMAP team's result (jKomats u et alJ 
120031 ISpergel et alJ 12006]) . In particular, the origin of 
the asymmetry in statistical properties of the CMB be- 
tween the_Gala£tic_jiOT^eni_^^ southern hemispheres 
l)Parkll2004t lEriksen et alJl2004bft has not been fully ex- 
plained. iFreeman et alJ l)2006j) have studied the effects 
of the map-making algorithm on the observed asymme- 
try in CMB temperatures between the Galactic north- 
ern and southern hemispheres in the WMAP data. An- 
other interesting issue is the low CMB quadrupole power. 
Cut-sky analysis of WMAP three-year data gives a 
quadrupole power of 211 /iK 2 , which is quite a bit lower 
than the ex pected value in the best-fit ACDM mode l 
(1250 /iK 2 ; IHinshaw et all 12001 ISpergel et all l200fi|) . 
Through subsequent analyses of the WMAP data, such 
a low quadrupole power has been confirmed and the ef- 
fect of the residual foreground emission on the statistical 
properties of CMB h a s been extensively studied (e.g., 
Te 



ies or (Jivia nas been extensively studied (e.g., 
■k et all 12001 lEfstathiol l200l lEriksen et alJ 



_ IBielewicx et al.ll20q4L boM ISlosar fc Seliald 
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2004 



200£ 
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Ide Oliveira-Costa fc Teginar k 2006). The alignment and 
planarity of I = 2 and 3 modes on the sky are also in- 
terest ing features that have been seen in the WMAP 
data ide Oliveira-Costa et aUl2004 iSchwarz et al.ll2004t 
|Lanci fc Magiieiioll200ft iConi et alJl200ft lAbramo et all 

Although each of the anomalies in the CMB anisotropy 
may have its own cosmological origin, there is a possibil- 
ity that the residual Galactic foreground emission has 
affected the nature of the observed temperature fluctu- 
ations. Therefore, foreground subtraction is very im- 
portant as a starting point of all CMB-related analy- 
ses. The most popular method of foreground removal 
is to model the Galactic emission as the weighted sum 
of foreground templates suc h as synchrotron, fre e-free, 
and dust emission maps (e.g. JBennett et al.l2003b(K An- 
other method of foreground removal is the internal lin- 
ear c ombination (LLC) method using multi-frequency 
data iBrandt et all 119941: iTegmark fc Efstathioul 11999: 



IBennett et al.ll2003bjh 

The WMAP team has created the ILC map (hereafter 
WILC1YR) by computing a weighted combination of the 
WMAP maps that have been band averaged within each 
of the five WMAP frequency ban ds, all smoothed to 1° 
resolution ijBennett et al.ll2003b|) . The weights are de- 
termined by minimizing the variance of temperatures in 
the combined map, with a constraint that the sum of the 
weights is one. They have defined twelve disjoint regions 
on the sky, within which weights are determined inde- 
pendently. Except for the biggest Kp2 region, the other 
11 regions are defined by subdividing the inner Galactic 
plane with lines of constant Galactic longitudes. How- 
ever, the combination weights found by the WMAP team 
do not give the minimum variance, and their method of 
'non-linear searching' for the set of weights has not been 
d escribed. 

lEriksen et all J2£ 



Mai) have improved the WMAP 
team's ILC method by applying a Lagrange multiplier 
to produce a variance-minimized ILC map (hereafter 
LILC). They use the same 12 disjoint regions and the 
smoothin g scheme to remove di s continuous bou ndaries 
as used in lBennett et alJ l)2003b|) . lEriksen et alJ has em- 
phasized that the effects of noise on the performance of 
the ILC method is very important. If noise is high, the 
ILC method finds the best combination weights that min- 
i mize the instrumental noise rather than foregrounds. 

ITegmark et al.l l|2003ft have applied a variant of the ILC 
method to make a foreground-cleaned CMB map (here- 
after TCM) with high resolution. By noticing that the 
Galactic foregrounds appear dominantly on larger an- 
gular scales while the instrument noise dominates only 
at smaller scales, they have applied linear combina- 
tions in harmonic space to remove foreground emission 
at different angular scales separately. Based on lev- 
els of the Galactic foreground intensity, ITegmark et alJ 
(2003) have defined nine disjoint regions, where fore- 
ground removal has been done indep endently. Recently, 
Ide Oliveira-Costa fc Tegmarkl (|200(tT> have made a new 
foreground-cleaned CMB map (hereafter TCM3YR) by 
applying the same technique to the WMAP three-year 
data. 

In the three- year data analysis, the WMAP team 
has made a new ILC map (hereafter WILC3YR) by 
combining the WMAP data at five different bands 



(Hinshaw et alJl200a) . From the one- year version of re- 
gion definition map, they have eliminated the Taurus A 
region that is too small to give a reliable foreground re- 
moval while they have added a new region to minimize 
the dust residuals in the Galactic plane (denoted as Re- 
gion 1 in Fig. Et). The y have used th e Lagra nge multi- 
plier method as used in lEriksen et alJ l)2004a[) to deter- 
mine the ILC coefficients. To quantify the bias due to 
the residual foreground emissions, they have performed 
one hundred Monte Carlo simulations using the one-year 
Galaxy model based on the Maximum Entropy Method 
(MEM). The three-year ILC map, WILC3YR, has been 
produced by subtracting this bias prediction. 

All ILC methods require several disjoint regions that 
contain map pixels with similar properties. Otherwise, 
the efficiency of the ILC method becomes very low, leav- 
ing significant residual foregrounds. The Galactic fore- 
ground emission appears as the sum of the individual 
emissions with different physical natures, and its spec- 
tral behavior varies over a wide range of frequencies and 
over the whole sk y. The disjoint region s defined by the 
WMAP team and lTegmark et alJ lj2003j) do not fully re- 
flect the Galactic foreground properties. 

In this paper, we derive a new foreground-reduced 
CMB map by applying the minimum-variance ILC 
method including information about the foreground 
spectral properties. The outline of this paper is as fol- 
lows. In §2, we describe our simple minimum- variance 
ILC method, and present a new definition of sky regions 
for ILC. In §3, we derive a new foreground-cleaned CMB 
map from the WMAP three-year data. Measurements 
of statistics that are related to the quadrupole and the 
octopole are given in §4. We discuss our results in §5. 

Throughout this paper, the power spectrum is calcu- 
lated by 

5T? = 1(1 + l)Ci/27r, (1) 



with 



Ci = (21 + 1)" 



E 



law. 



(2) 



where a/ m are the spherical harmonic coefficients in the 
Galactic coordinate system. 

2. INTERNAL LINEAR COMBINATION METHOD 

2.1. ILC with a Minimum Variance Constraint 

We find an optimal weight set for linear combination of 
observed CMB maps to remove the Galactic foregrounds 
with the constraint that the combined map have the min- 
imum variance. Let us denote the observed temperature 
at p-th pixel and at i-th frequency band by Ti(p). The 
observed temperature is a sum of the CMB tempera- 
ture, instrument noise, and the Galactic emission. We 
perform ILC on a region 1Z by minimizing the variance 
of the weighted sum of observed temperatures at pixels 
within the region. For multi-frequency CMB data with 
k bands, the temperature of p-th pixel of the linearly 
combined map can be written as 

k fc-1 
T(p) = ]T wtTiip) = WiD ik {p) + T k (p), (3) 



i=i 



= 1. The 



where D ik (jp) = Tj(p) - T k (p) and £\ = 
variance of T(j>) over pixels within the region 1Z is given 




Fig. 1. — The three-year version of Galactic foreground emission maps at five WMAP frequencies (K, Ka, Q, V, and W bands) made 
by co-adding synchrotron, free-free, and dust emission components derived from the MEM. The resolution of each map is 1° FWHM. 
The Mollweide projection in Galactic coordinates is used to display all maps. In each figure, the Galactic center is located at the center 
and the Galactic plane lies at the horizontal line that passes through the Galactic center. Galactic longitude increases from center to left 
(I = 0°-180°) and from right to center (I = 180°-360°) while Galactic latitude increases from bottom to top (6 = -90°-+90°). 



by 

a 2 (T) 



pen 
fc-i 



pen 



-, 2 



J2 W {Afc(p) - Afc} + {Tk(p) - fk) 



(4) 



where T represents the average of T(p) over Nn pix- 
els within the specified region 1Z. ILC coefficients that 
give the minimum variance is obtained by solving linear 
equations da 2 (T) / dwj — 0, where j = 1, - • • ,k — 1. The 
solution w = (w±, • • • , Wk-i) T is simply obtained from 

w = -S -1 t. (5) 

Here components of the (k— 1) x (k — 1) symmetric matrix 
S and the (k — 1) x 1 vector t are given by 
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where i, j = 1, ■ • • , k — 1. The weight for the last fre- 
quency band map is Wk = 1 — J^^i w i- With the set 
of ILC weights, a foreground-reduced CMB map can be 
obtained from equation J2J (see SQ. 

2.2. Defining Pixel-Groups with Common Foreground 
Spectral Properties 

The efficiency of foreground reduction by the linear 
combination method becomes highest when the region 
contains pixels with the same foreground spectral varia- 
tion over frequencies. On the other hand, if a region con- 
tains a large number of pixels whose Galactic foreground 
spectral property varies dramatically over the sky, as is 



the region K = of WILC3YR (see Fig. gt), the ILC 
just tries to find the best-fit combination weights that is 
optimal only to the average foreground spectral behav- 
ior. The ILC then results in positive or negative residual 
bias in the map (3S3J- 

To assess the Galactic foreground properties, we use 
the MEM Galactic foreground maps derived by the 
WMAP team, which are shown in Figure^ The WMAP 
team has tried to distinguish different emission sources 
from one another by applying the MEM to the WMAP 
data, where the prior spatial distribution and spectral 
behavior of foreground components have been assumed 
by using the Galactic template maps. From this informa- 
tion the team has produced synchrotron, free-free, and 
dust emission maps with 1° resolution at each WMAP 
frequency. For each band, we co-add three emission maps 
to make a Galaxy foreground map, upgrade its resolu- 
tion to Nside = 512, and further smooth the map with a 
Gaussian filter with variable widths (0° - 2°) to reduce 
the MEM reconstruction noise. We smooth each Galaxy 
map over FWHM = 0?5, and then subtract it from the 
unsmoothed Galaxy map to make a difference map, from 
which the standard deviation is calculated at the region 
with |6| > 60°. The signal-to-noise ratio at each pixel 
for each band is defined as the foreground intensity di- 
vided by the standard deviation. The minimum value 
among the five signal-to-noise ratios at each pixel is con- 
sidered as the representative signal-to-noise ratio (S/N). 
The width of the smoothing filter has been set to zero for 
S/N > 10, and FWHM = 0.2(10 - S/N) for S/N < 10. 

Tabled lists the spherical harmonic coefficients (a; m ) 
of the K-band Galactic foreground map up to / = 6. 
The aim's of foreground emissions at higher frequency 
bands have similar patterns with smaller amplitudes. 
The Galactic foreground emission is strong at even l- 
modes, among which I = 2 and m = mode is the 
strongest one. The real components of even (odd) m- 
modes for even (odd) I are stronger than the imaginary 
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TABLE 1 

low /-mode spherical harmonic coefficients of the 
K-band Galactic foreground emission 
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components. Modes of strong Galactic signal are ex- 
pected to affect the corresponding modes of CMB sig- 
nal significantly, resulting in biases in the reconstruction 
maps. 

Figure [2t shows intensity variations of the Galactic 
foreground emission across WMAP bands at about 200 
randomly selected pixels. Figure^ shows the histograms 
of foreground spectral indices (a) within K-Ka, Ka-Q, 
Q-V, and V-W frequency intervals measured at all pixels 
in the sky. Each spectral index is measured by modeling 
the foreground intensity by I v oc v a . The spectral indices 
measured at frequency interval V-W (dot-dashed curve) 
have a wide range of values from negative to positive, 
centered at a w —0.6 unlike those at K-Ka, Ka-Q, and 
Q-V that are located at a < — 1 (dotted, dashed, and 
long-dashed curves). Therefore, it is reasonable to treat 
spectral indices at low (K-V) and high (V-W) frequency 
intervals separately. For simplicity, we use an average 
spectral index measured for four foreground intensities 
from K to V bands as the representative spectral index 
at low frequency bands (Fig. solid curve). Hereafter 
we call the spectral indices measured at the low and high 
frequency intervals a and /3, respectively. 

Distributions of the spectral indices on the sky mea- 
sured at low and high frequency intervals are shown in 
Figure |2| The Galactic plane region has spectral indices 

of a ps -2.5 2.0 and /3 « -1.0 - +1.0 while the 

high latitude region with low foreground contamination 
shows large fluctuations, especially for (3 map. ft should 
be noted that the spectral index is a sensitive function 
of both frequency and direction, and that the foreground 
subtraction of 1LC will be most successful when the ILC 
is performed at a region where the foreground character- 
istic is most homogeneous. 

For each of the two histograms of spectral indices at the 
high and low frequency bands (Fig. [20, we define twenty 
spectral index bins so that each bin contains equal num- 
ber (157, 300) of pixels with similar foreground spectral 
properties. By combining the two sets of spectral-index 
bins, one at low and the other at high frequency band, 
we define four hundred (20 x 20 = 400) groups of pixels 
with similar spectral indices. Because pixels within the 
same a-bin have a wide range of spectral index /3, the 
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Fig. 2. — (a) Intensity variation of the Galaxy foregrounds over 
WMAP frequencies (22.8 [K], 33.0 [Ka], 40.7 [Q], 60.8 [V], and 
93.5 GHz [W]). According to levels of Galactic latitudes (low to 
high), plots are shown with different colors (red to blue). Only 
about 200 pixel values are shown to avoid overlaps, (b) Histograms 
of foreground spectral indices measured at K-Ka (dotted), Ka-Q 
(dashed), Q-V (long-dashed), and V— W (dot-dashed curve) fre- 
quency intervals. The solid curve denotes a histogram of spectral 
indices measured from four foreground intensities at K to V bands. 
The MEM Galactic foreground maps smoothed with a Gaussian 
filter of variable widths have been used for spectral index measure- 
ment. 



number of pixels contained in each group differs from 
group to group. The resulting group index map is shown 
in Figure 0Jt. The variation of the spectral index across 
the sky is naturally taken into account in this way. We 
note that pixels at the high latitude regions are assigned 
with a wide range of group indices. The ILC will be ap- 
plied separately for each group of pixels with similar fore- 
ground spectral properties to obtain foreground-reduced 
CMB temperatures in the next section. 

3. APPLICATION TO THE WMAP THREE- YEAR DATA 

3.1. The ILC Map Derived From the WMAP 
Three- Year Maps: SILC400 
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Fig. 3. — All-sky maps showing distributions of foreground spec- 
tral indices measured at (a) low (K-V) and (b) high (V-W) fre- 
quency intervals. The MEM Galactic foreground maps smoothed 
with a Gaussian filter of variable widths have been used to measure 
the spectral indices. 



We apply our ILC method to the WMAP three-year 
data 1 . The WMAP has one K band (22.8 GHz), one Ka 
band (33.0 GHz), two Q band (40.7 GHz), two V band 
(60.8 GHz), and four W band (93.5 GHz) differencing 
assemblies, with 0?82, 0?62, 0?49, 0?33, and 0?21 FWHM 
beam widths, respectively. The WMAP m aps are made 
in the HEALPix 2 format with Nside = 512 ijGorski et all 
119991 12005(1 . The total number of pixels of a map is 
12 x Nside 2 = 3, 145, 728. 

We use the five band-averaged WMAP maps with 1° 
FWHM resolution that are presented by the WMAP 
team. The maps have been produced in the following 
way. Each map of each differencing assembly is decon- 
volved with the corresponding channel-specific beam and 
convolved with FWHM = 1° beam. The maps of the 
same frequency band are averaged with noise-weight at 
each pixel taken into account, and finally the five WMAP 
maps at K, Ka, Q, V, and W bands with FWHM = 1° 
resolution are obtained. 

We obtain a foreground-reduced ILC map (hereafter 
WILC12) from the three-year WMAP data with 1° 
FWHM resolution at five frequency bands (k — 5) by 
using equations ©, © and © for each of 12 disjoint 
regions that are defined by the WMAP team (Figs. Efr 
and b). Our weig ht coefficients are sli ghtly different 
from those found bv lHinshaw et all l|2006|) who compute 
minimum-variance ILC weights by using Lagrange mul- 

1 http://lambda.gsfc.nasa.gov 

2 http://hcalpix.jpl.nasa.gov 



tipliers, but our weight sets give smaller variances. To 
estimate the residual bias in the WILC12 map we gen- 
erate two hundred maps mimicking the WILC12 map by 
combining theoretical CMB temperature signal with in- 
strument noise and foreground emissions. The three-year 
version of MEM Galactic foreground maps have been 
used as the contaminating sources. The average and the 
standard deviation maps of residual emission in two hun- 
dred WILC12 simulation maps with respect to the true 
CMB temperatures are shown in Figures^: and d. From 
the WILC12 simulations, we found that ILC applied to 
the large high latitude region that contains pixels with a 
wide range of foreground spectral indices (e.g., the region 
with index 1Z = 0) induces residual biases whose levels 
are different at different positions on the sky. As shown 
in Figure Et, the residual emission is negative near the 
Galactic plane and is positive at high latitude region. 

Figure 0t shows a new foreground-reduced CMB map 
obtained by applying the ILC method described in i )2.1l 
with the minimum variance constraint, where the com- 
binations are performed for four hundred pixel-groups 
of common spectral properties independently (hereafter 
SILC400). The SILC400 shown here is a map that has 
been corrected for the residual bias based on the ILC 
simulation as described in H3.2I and the monopole and 
dipole have been subtracted from the map. The map 
has been shown after further smoothed by a Gaussian 
filter of FWHM = 1° to reduce discontinuities between 
boundaries of pixel groups. In the subsequent analysis 
of large scale modes, however, we use the unsmoothed 
SILC400. 

ILC combination weights for 400 groups are plotted in 
Figure El The combination weights varies dramatically 
from group to group, except for the K-band weights (wi; 
filled circles) that are stably distributed around —0.25. 
Most of the weights in Ka-W bands (1V2-W5) are positive 
for group index < 250 corresponding to high latitude 
regions, while they are fluctuating at the Galactic plane 
regions (group index > 250). 

Comparisons of SILC400 with other foreground- 
reduced maps (WILC3YR and TCM3YR) are shown as 
difference maps in Figure where all maps are made 
to have 1?414 FWHM resolution with the monopole and 
the dipole components removed before differencing. The 
WILC3YR is similar to SILC400 within 40 difference, 
except for the Galaxy center and cirrus regions. The 
cirrus is apparent as red and the Galactic center looks 
partly blue. For TCM3YR, its difference from SILC400 
shows the region boundaries near the Galactic plane and 
mid-latitude, especially at I — 0° - 180° (left part in Fig. 
Ep) , where positive and negative temperatures appear in- 
terchangeably. Those region boundaries were defined by 
iTeemarket aDl|200ll) . 

Low /-mode spherical harmonic coefficients (a/ m ) of 
SILC400 are listed in Table (I = 2 coefficients for 
WILC3YR and TCM3YR are summarized in the note). 
We note that the SILC400 has a 20 (= 7.51 /jK) that 
is between those of TCM3YR (3.34 fjK) and WILC3YR 
(11.5 /iK), and also has slightly different 022- The ampli- 
tudes of 020 and 022 are bigger than 021 both in SILC400 
and the Galactic foreground map (Tabled . However, 020 
and 022 have different signs in two maps: the quadrupole 
of CMB and the Galaxy are negatively correlated (see 
|Q. This is mainly due to the big cold spot in the CMB 
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Fig. 4. — (a) A group index map showing four hundred groups of pixels with similar foreground properties denned by using the foreground 
spectral index maps at low and high frequency bands (Fig. The group index runs from to 399. (ft) The foreground-reduced CMB 
temperature fluctuation map (SILC400) derived from the five-band WMAP maps by applying the ILC method which adopts the pixel- 
groups as defined in (a). The SILC400 has been corrected for the residual bias, and the monopole and dipole have been removed from the 
map. Due to the further smoothing with a Gaussian filter of FWHM = 1°, the map resolution is 1?414 FWHM. The SILC400 is available 
at http://newton.kias.re.kr/~parkc/CMB/SILC400.html together with its spherical harmonic coefficients ai m up to Z = 200. 



map at the Galactic center region (see Fig. . 

3.2. Simulations: Test for the ILC Method 

As shown in SILC400 differs from foreground- 

reduced CMB maps previously made by others. Because 
the ILC method conserves the blackbody nature of the 
CMB signal, all maps contain exactly the same informa- 
tion for CMB temperature fluctuations but with different 
level of foreground residuals. To quantify the level of the 
residual foreground in our variance-minimized ILC map, 



we have performed two hundred simulations that mimic 
the WMAP data and analyzed them in the same way 
that SILC400 is made. 

First, we simulate two hundred WMAP Gaussian 
CMB signal maps for each differencing assembly at each 
band. The concordance flat ACDM model power spec- 
trum that fits to th e WMAP data only, has been used 
ijSpergel et al.ll2006|) . For each mock observation, the 
same aim's are used for all frequency channels. During 
the map generation, the WMAP beam transfer function 
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Fig. 5. — (a) Region definition map defined by the WMAP team (12 regions), (b) A foreground-reduced CMB map derived with the 
minimum-variance ILC method (WILC12). The ILC has been applied to 12 regions independently. Maps showing (c) the average and (d) 
the standard deviation of the residual foreground emissions in the ILC map estimated from two hundred WILC12 simulation maps. The 
monopole and dipole components of the average map have been removed. 
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Fig. 6. — Minimum- variance ILC weights of SILC400 obtained 
for four hundred pixel-groups. The ILC weights are denoted as 
yellow filled circles for K (wi), red open circles for Ka (W2), green 
open squares for Q (ws), blue filled squares for V (104), and violet 
crosses for W (1^4) bands. 



Bi is used for each differencing assembly, and the instru- 
ment noise at each pixel is randomly drawn from the 
Gaussian distribution with variance of o^/Nobs, where 
7V b s is the effective number of observations at each pixel, 
and ctq is the global noise level of the map (Table 1 of 



TCM3YR - SILC400 

Fig. 7. — Comparison of SILC400 with other foreground-reduced 
CMB maps (a) WILC3YR and (b) TCM3YR. Shown are differences 
of WILC3YR and TCM3YR from SILC400. Before differencing, 
all maps are prepared to have FWHM = 1?414 resolution and the 
monopole and dipole modes are removed. 



TABLE 2 

LOW /-MODE SPHERICAL HARMONIC COEFFICIENTS OF THE 

Minimum- Variance ILC Map SILC400 



I 


m 


■ — 

Re[a; m ] 




Im[a !m ] 


I 


in 


■ — 

Re[a Im ] 




Im[a im ] 


2 





7.51 


0.00 


5 





14.33 


0.00 




1 


-1.54 


4.82 




1 


23.96 


1.09 




2 


-18.54 


-18.01 




2 


-7.75 


2.78 


3 





-5.25 


0.00 




3 


19.75 


3.58 




1 


-10.74 


3.76 




1 


-4.21 


8.53 




2 


22.54 


2.13 




5 


9.48 


20.26 




3 


-13.52 


30.10 


6 





-1.76 


0.00 


4 





23.31 


0.00 




1 


0.78 


4.17 




1 


-6.97 


9.12 




2 


10.39 


-2.75 




2 


11.14 


4.26 




3 


-9.02 


-0.28 




3 


7.34 


-21.02 




1 


13.73 


-1.56 




1 


1.75 


-9.68 




5 
6 


-5.07 
3.64 


-3.74 
8.91 



Note. — The coefficients are calculated in Galactic 
coordinates in units of fiK. The / = 2 coefficients for 
WILC3YR are a 20 = 11.48, o 2 i = -0.05 + 4.86i, and 
a 22 = -14.41 - 18.80i. For TCM3YR, a 20 = 3.34, 
a,2i = 0.26 + 4.88i, and a 22 = -14.88 - 17.26i. 



iHinshaw et aHl2006|) . 

Secondly, for each differencing assembly realization, 
we deconvolve the beam effect, and convolve the map 
to a common resolution of 1° Gaussian FWHM. Differ- 
encing assembly maps at the same band are averaged 
with weights given by the noise variance at each pixel 
to produce an average map at each frequency band. Fi- 
nally, the Galactic foreground map (with 1° FWHM res- 
olution) at each frequency is co-added to the average 
WMAP CMB map. We do not use the Galactic tem- 
plate meth od to mimic th e Galact ic foreground emissions 
as used in lEriksen et alJ l)2004a|) . but use the Galactic 
foreground maps derived by MEM from the three-year 
WMAP data. Compared to weighted average of template 
maps, these foreground maps reasonably describe the 
Galactic foreground emission even at the Galactic plane. 
However, the Galactic foreground maps derived by the 
MEM are somewhat noisy because of the MEM recon- 
struction noise. For each simulation, the five foreground- 
added CMB maps are put into the same ILC pipeline as 
used in A3. II 

Figure [S] shows average and standard deviation maps 
of the residual foreground emissions with respect to true 
CMB temperatures in two hundred SILC400 simulation 
maps. It demonstrates that except for the Galactic 
plane region our ILC method correctly reconstructs the 
CMB signal without significant bias due to the Galactic 
foreground contamination unlike the average residual- 
emission map of WILC12 simulations (Fig. How- 
ever, some residual foreground features with 10-30 /xK 
levels are seen at the high Galactic latitude regions, 
where the foreground emission is very weak at high fre- 
quency bands (V and W) and the spectral index (3 has 
been measured with large uncertainty due to the MEM 
reconstruction noises (upper right part of Fig. IBfr). The 
standard deviation map indicates that SILC400 may be 
contaminated by the residual foreground at the Galactic 
plane and the mid-latitude regions. Compared with the 
WILC12 case (Fig. [SJQ, the standard deviation map has 
larger amplitudes at high latitude region, because pixel- 





Fig. 8. — (a) Average and (6) standard deviation maps of the 
residual foreground emissions with respect to true CMB tempera- 
tures in two hundred SILC400 simulation maps. 



groups used for SILC400 contain much fewer pixels than 
the TZ = region of WILC12. 

Figures Eh, 6, and c compare the spherical harmonic 
coefficients of I — 2 mode calculated from two hundred 
SILC400 simulation maps with those of true input CMB 
maps. To reduce the effect of the residual bias, we have 
subtracted the average map of the residual foreground 
emission from each SILC400 simulation map before cal- 
culating the ai m ■ The (120 and the real component of 022 
are contaminated by the residual Galactic emission in 
the sense that each plot, when fitted with a straight line 
(dashed lines), gives a slope less than 1, and some scatter 
is seen in the 020 plot. Quadrupole powers from two hun- 
dred SILC400 simulations against true CMB quadrupole 
powers are shown as open circles in Figure^, where the 
quadrupole power predicted in the concordance ACDM 
model is denoted a s horizontal line (5Tf = 1250 (1K 2 ; 
ISpergel et all [20061. The SILC400 systematically un- 
derestimates the CMB quadrupole powers, but with 
much less scatter th an LILC simulations (see Fig. 6 of 
Erik sen et aT]l2004a|) . 

Spherical harmonic coefficients for I — 3 and 4 of 
SILC400 simulations are compared with true values in 
Figure^D The I — 3 modes are excellently reconstructed 
in the SILC400 simulations, although plots for real com- 
ponents of 0,31 and 0,33 against true values show slight 
decrease in correlation slopes. For I — 4, the real compo- 
nents of 040 and 042 show similar patterns to I — 2 case. 
The residual foreground emission has induced slope de- 
crease and scatter in the 040 plot (Fig. ITTJfe). 

We can remove the systematic effects of the ILC fore- 
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a„ (**K, input CMB) (,/j.K. input CMB) 




(,uK, input CMB) <3T= (juK 3 , input CMB) 



Fig. 9. — (a)-(c) Plots of / = 2 spherical harmonic coeffi- 
cients (aj m ) calculated from two hundred SILC400 simulation maps 
against those from true CMB maps (filled circles for real and the 
open circles for imaginary components). For real components of 
even m modes, each plot has been fitted with a straight line (dashed 
lines), (d) Plot of quadrupole powers measured from SILC400 sim- 
ulation maps against those from the true CMB maps (open cir- 
cles). Quadrupole powers measured from bias-corrected a; m 's of 
SILC400 simulations maps are plotted against the true values as 
filled circles. 



ground reduction method on ai rn of the SILC400 simula- 
tion maps. We have estimated the true ai m up to I = 200 
by applying a simple relation a; m = (a} m c —b)/a, where 
a }m ^ s the spherical harmonic coefficient obtained from 
each SILC400 simulation map from which the average 
residual bias map of Figure |Bb has been subtracted, and 
a and b denote a correlation slope and an offset from the 
line-fit. Linear fit parameters have been derived from the 
aim plots, where each slope and offset are obtained for 
real or imaginary components of aim's, independently. 
Real components of ai m for even (odd) I and even (odd) 
to modes show the systematic decreases in the slope of 
the linear relation. However, in all cases, no large offset is 
observed in the SILC400 simulations. The result for I = 2 
is shown as filled circles in Figure El^- The quadrupole 
powers of SILC400 simulation maps, which were under- 
estimated before correction, now have the correct mean 
amplitudes compared with the true values. Figure If fl 
shows the low i-mode powers of SILC400 simulations ver- 
sus true CMB powers for I = 3-6 before (open circles) 
and after (filled circles) the bias correction for ai m . 

4. STATISTICS OF LARGE-SCALE MODES 

Table summarizes the quadrupole and octopole pow- 
ers measured from the WMAP data in the previous stud- 
ies, including our new estimates. For S1LC400 results, 
the 68% confidence limits have been deduced from the 
distribution of <5T 2 (ILC) - ST 2 (CMB) in two hundred 
SILC400 simulations. The third column of Tabic lists 
the probabilities (j>) of the quadrupole power being as 



low as me asured value if the concordanc e ACDM model 
is correct (|de Oliveira-Costa et al.ll2~004ft . However, the 
p values for S1LC400 maps have been estimated from 

p = J dQ ohs p(Q < QobsMQobs), (7) 

where Q denotes quadrupole power and p{Q < Qobs) is 
the probability of finding the quadrupole power as low 
as the observed value Qobs- The p(Q bs) is the prob- 
ability of observing the quadrupole power Q bs and is 
calculated from the distribution of quadrupole powers 
from two hundred SILC400 simulations. The quadrupole 
power of SILC400 (244.5±^ 3 8 3 fiK 2 ) is close to that 
of WILC3YR (248.5 fiK 2 ) and is between values of 
TCM3YR (209.6 fiK 2 ) and LILC (350.6 pK 2 ) . Our esti- 
mates compare with those of lBielewicz et ahl ll2004|) who 
have applied power equalization filter to the high lati- 
tude WMAP data and obtained <5T 2 2 = 279 ± 49 and 
ST 2 = 870 ± 98 from TCM. 

Also listed in Table |3| are the quadrupole and oc- 
topole powers from S1LC400 after removing the ex- 
pected biases in ai m due to the minimum variance 
ILC method. Our bias-corrected quadrupole power 
(275.8± ^ n nK 2 ) is higher than the WMAP team's mea- 
surement ( Hinshaw et al. 2006) , but they are still consis- 
tent with each other. The probability of observing such 
a low quadrupole power in the concordance ACDM uni- 
verse is as high as p — 5.7%. 

Figure ^] compares S1LC400 power spectrum with 
other measurements up to / = 35 (top) and I — 200 (bot- 
tom). Also shown in the top panel is the power spec- 
trum of K-band MEM Galactic foreground map whose 
zigzag pattern at I = 2-8 is very similar to that of the 
CMB power spectra. For WILC3YR, TCM3YR, and 
SILC400, we have deconvolved the beam effects to ob- 
tain correct amplitudes of power spectrum. The WMAP 
team's power spectrum measurement from the cut-sky 
analysis and the best-fit ACDM power spectrum are also 
shown for comparison with 68% and 95% con fidence 
limits l|Hinshaw et all 120061 iSpergel et all 12006ft . The 
quadrupole power of the ACDM model is 1250 /iK 2 with 
68% and 95% confidence intervals, 510—1961 /iK 2 and 
198-3311 /zK 2 , respectively (fr om the ten thousand real- 
izations) . 

Our new estimate of power spectrum from S1LC400 
(red curve) is very similar to others up to I « 30. Due to 
the instrument and reconstruction noises, all ILC maps 
give angular power spectra with higher amplitudes with 
increasing Z, compared with the ACDM model. As shown 
in Figure^] (bottom), the angular power spectrum of the 
SILC400 has lower amplitude than those of W1LC3YR 
and TCM3YR up to I = 200, which indicates that our 
ILC method better satisfies the minimum variance con- 
straint. 

It has been reported that both the quadrupole and 
the octopole of the WMAP map appear planar, with 
most of their hot and cold spots located on a single plane 
in the sky, and the two plan es appear roughly aligned 
l)de Oliveira-Costa et alJ20T)4T) . A simple way to quantify 
a preferred axis for arbitrary multipoles is to find the axis 
n around which the angular momentum dispersion 

^|(n.L) 2 |^)=£m 2 K m (n)| 2 (8) 
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Fig. 10. — Plots of spherical harmonic coefficients (a; m 's) estimated from two hundred S1LC400 simulation maps against those from true 
CMB maps for (a)— (d) I = 3 and (e)— (ft.) I = 4 modes. The 043 plot is omitted because it is very similar to 041 plot. 



TABLE 3 

CMB QUADRUPOLE AND OCTOPOLE POWERS 



Measurement 


5T| (p.K 2 ) 


p value 


«5T| (/1K 2 ) 


Spergel et al. best-fit ACDM model 


1250 




1143 


Hinshaw et al. cut sky analysis 


211.0 


2.6% 


1041 


WMAP team's ILC map (WILC3YR) 


248.5 


3.7% 


1051 


Tegmark et al. cleaned map (TCM) 


201.6 


2.3% 


866.1 


de Oliveira-Costa & Tegmark cleaned map (TCM3YR) 


209.6 


2.5% 


1037.8 


Eriksen et al. ILC map (LILC) 


350.6 


7.6% 


1090 


Minimum-variance ILC map (SILC400) 


244 e;+ 8 3- 8 


4.9%t 


859A -tofo 


SILC400 (bias corrected) 


275 R+94-3 


5.7%t 


951.9_ g3 Q 



Note. — p value denotes the probability of finding a lower qua drupole than the 
measured value if the concordance ACDM model (Socracl ct al. 200(|) is correct. The p 
values with f have been estimated from equation Q . 



is maximized i|de Oliveira-Costa et al.ir2004|) . where the 
CMB map is considered as a wave function AT(n) = 
ip(n). Here ai m (h) denotes the spherical harmonic coef- 
ficients of the CMB map in a rotated coordinate system 
with its z-axis in the n direction. For each I we find a unit 
vector (fii) that maximizes the angular momentum dis- 
persion. We start to evaluate equation © for all the unit 
vectors corresponding to HEALPix pixel centers at res- 
olution Nside = 32, perform the same operation only for 
pixels at higher resolution around the unit vector found 
in the previous step, and finally obtain a unit vector that 
maximizes equation (|SJ) at resolution Nside = 4096. 

Table 0] lists directions of the quadrupole and the 
octopole and separations between two poles measured 
from the WMAP data. The 68% confidence limits 
for the measured angular separation 6*23 of SILC400 
have been estimated from the distribution of 6*23 (ILC) — 
#23 (CMB) in two hundred SILC400 simulations. Since 
the quantity \fi2 • 1^3 1 has a uniform distribution on 
the unit interval [0,1], the probability of finding an 
angular separation smaller than the measured separa- 
tion is simply given by p = 1 — \fi2 ■ T1.3 1 (last col- 



umn in Table EJ. However, the p values for SILC400 
maps have been calculated from equation J7J with Q 
replaced with #23- The bias-corrected SILC400 has an 

angular separation (#23 = 11?8 + i*J that is somewhat 

larger than those from TCM, LILC, and WILC3YR 
but is very similar to that of TCM3YR. We find that 
the octopole directions (1^3) are more stable than the 
quadrupole directions (1^2) ■ Our result for SILC400 
confirms the previous results tha t the quadrupo l e and 
octopole dire ctions a re aligned iTegmark et all l2003t 
[ Eriksen et all l2004at fBielewicz et all 120041 iCopi et alJ 
l2004ll2006HSchwarz et al.ll2004[) . with p = 4.3%. 

Figure 1131 compares angular separations between 
quadrupole and octopole directions of two hundred 
SILC400 simulation maps with the true values, and 
demonstrates that SILC400 can reconstruct the correct 
multipole directions of low /-modes. 

We also measure the degree of planarity of low 
I modes by calculating the t statistic defined by 
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TABLE 4 

Directions of the Quadrupole and the Octopole 



Maps 


h 


b 2 


h 


b S 


623 


p value 


WMAP team's ILC map (WILC3YR) 
Tegmark et al. clean map (TCM) 
TCM3YR (de Oliveira-Costa & Tegmark 2006) 
Eriksen et al. ILC map (LILC) 


235? 23 
257? 58 
224? 04 
247? 52 


68? 42 
58? 82 
76? 56 
61? 90 


236? 94 
238? 37 
236? 75 
232? 84 


62?61 
62? 04 
63?98 
63?41 


5? 85 
9? 96 
13? 23 
6? 89 


0.52% 
1.51% 
2.65% 
0.72% 


Minimum-variance ILC map (SILC400) 


240? 78 


75? 24 


234? 47 


62? 26 


id - ib _ 6 ?79 


5.36%+ 


SILC400 (bias corrected) 


245? 54 


74? 38 


234? 58 


63? 25 


11 077+6? 40 
lX - -8? 03 


4.29%t 



Note. — 623 is the angular separation defined as 623 = cos 1 (fi2 ■ ns), where n 2 — (£2,62) and 
1^3 = (£3,63). The p values with f have been estimated from equation Q with Q replaced by #23. 




1000 3000 1000 2000 3000 

6Tf input CMB) fflj {fiV?, input CMB) 

Fig. 11. — (a)—(d) Plots of low i-mode powers (/ = 3—6) estimated 
from two hundred SILC400 simulation maps against those from 
true CMB maps before (open circles) and after (filled circles) the 
effect of bias due to the residual Galactic emission on a; m 's is 
removed. 
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The t statistic measures the maximal percentage of l- 
mode power that contributes to |m| = I, and we obtain t 
value for each I by finding a direction n that maximizes 
equation Q. The maximization is performed over pix- 
els of a map in the similar way as used in finding the 
quadrupole and octopole directions. The performance of 
our ILC method in reconstructing the true t-statistic is 
shown in Figure^! where t values of SILC400 simulation 
maps against those from true CMB maps are plotted for 
/ = 2-6 modes. 

The t values measured from the foreground-cleaned 
maps are summarized in Table |SJ Numbers in the paren- 
theses represent the number of occurrences that have 
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Fig. 12. — Power spectra measured from SILC400 (thick red), 
WILC3YR (blue), and TCM3YR (green curve). The power spec- 
trum of K-band MEM Galactic foreground map has also been 
shown for comparison (brown curve) . Beam effects are deconvolved 
for the correct amplit udes of power spectra . The WMAP power 
spectrum measured by Hinshaw et al. 1200(f) (grey curve), and the 
ACDM model of Socrecl et al. 12006) (dashed curve) are shown 
for comparison with 68% and 95% confidence limits (grey dashed 
curves) . 



larger t values than the measured value among the two 
hundred SILC400 simulations. The 68% confidence lim- 
its for I = 2 and 3 are obtained from the distributions of 
i(ILC)-i(CMB) in the SILC400 simulations. The t value 
for I = 3 mode of SILC400 (bias corrected; t = O^lt^) 
is very similar to the previously measured values, and im- 
plies that the octopole has high planarity. However, the 
planarity is not statistically significant because the prob- 
ability that SILC400 simulation map has t value larger 
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Fig. 13. — Plots of angular separation (623) between the 
quadrupole and the octopole directions measured from two hun- 
dred SILC400 simulation maps against those from true CMB maps, 
before (open circles) and after (filled circles) the bias correction. 
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Fig. 14. — Plots of i-statistic for I = 2—6 modes measured from 
two hundred SILC400 simulation maps against those from true 
CMB maps, before (open circles) and after (filled circles) the bias 
correction. 



than the measured t value is 18.5%: thirty seven among 
two hundred SILC400 simulations. The probability of 
having the t value for the quadrupole mode as high as 
the measured value of SILC400 (t = 0.98±£^) is still 
higher (98/200). As shown in Figure ITU most of the 
CMB quadrupole modes have t > 0.8. 

As pointed bv lEriksen et all l)2004af) . the I = 5 and 
I = 6 modes are very peculiar in their symmetry prop- 
erties. Only one case out of the two hundred SILC400 
simulations has t values larger than t = 0.38 for I = 5, 
and only four have t values larger than t = 0.83 for I = 6, 
which indicates that the distribution of temperature fluc- 
tuation of I = 5 mode is highly symmetric, and the I — 6 




1 10 100 1000 10' 1 10 100 1000 10' 

Number of Groups Number of Groups 

Fig. 15. — Standard deviations (cr; top-left), the quadrupole and 
the octopole powers (5T ; 2 ; bottom-left), and angular separations 
between the quadrupole and the octopole (#23; bottom-right panel) 
measured from ILC maps made for various number of pixel-groups 
from 4 to 6400. The top-right panel shows the variation of the 
quadrupole (filled circles) and the octopole (open circles) directions 
in Galactic coordinates. The horizontal line in the top-left panel 
denotes the standard deviation of the WILC3YR map. 



mode is planar together with the I = 2 and I = 3 modes. 

We also investigate whether the measured statistics de- 
pend on the number of pixel-groups and the smoothing 
scale. First, we measure standard deviation, quadrupole 
and octopole powers, and angular separation between 
the two multipoles from the ILC maps made by fixing 
the resolution to 1° FWHM while varying the number of 
pixel-groups in the group-index map from 4 (2 x 2) to 
6400 (80 x 80). The result is shown in Figure [U The 
standard deviation and quadrupole power decrease with 
increasing number of pixel-groups, while the octopole 
power is very stable. If the number of pixel-groups is 
larger than 100, all statistical quantities become stable. 
However, when the number of pixel-groups is small, a 
wide range of foreground spectral indices is allowed in 
a common group, and the ILC results in higher values 
of standard deviation, quadrupole power, and larger an- 
gular separation 023. We note that 6*23 becomes large 
because the quadrupole direction is unstable in this case 
(filled circles in the top- right panel of Fig. I15|) . 

Secondly, we measure the large scale mode statistics 
from the ILC maps produced from the WMAP data and 
the Galactic foreground maps with various angular res- 
olution. Each WMAP map or MEM-derived Galactic 
foreground map, which has originally 1° FWHM resolu- 
tion, has been further smoothed with a Gaussian filter of 
FWHM = 0?5-3?0 with steps of 0?5 (the total resolution 
varies from 1?12 to 3? 16 FWHM). The smoothed maps 
are used to calculate spectral indices for pixel-group defi- 
nition and for subsequent application to the ILC pipeline. 
In all cases, the number of pixel groups has been fixed 
to 400. As shown in Figure our ILC method gives 
stable and consistent amplitudes of the quadrupole and 
the octopole powers, although the angular separation 823 
tends to increase with increasing smoothing scale. 

5. DISCUSSION 

In this paper, we have derived a new foreground- 
reduced CMB map by applying a simple internal lin- 
ear combination method to the WMAP three-year 
data. Rather than using the disjoint sky regions as in 
WILC3YR, TCM, and LILC, we have defined a new 
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TABLE 5 

t Values of Low Spherical Harmonic Modes 



Maps 


I = 2 




/ = 3 




I = 4 


I = 5 


I = 6 


WMAP team's ILC map (WILC3YR) 
Tegmark et al. clean map (TCM) 
TCM3YR (de Oliveira-Costa & Tegmark 2006) 
Eriksen et al. ILC map (LILC) 


0.930 
0.957 
0.995 
0.913 




0.917 
0.942 
0.877 
0.934 




0.648 
0.588 
0.679 
0.587 


0.375 
0.372 
0.358 
0.374 


0.798 
0.783 
0.749 
0.806 


Minimum-variance ILC map (SILC400) 
SILC400 (bias corrected) 


u.yyi_ 

n q77 +0.016 
u - 9 "-0.024 


(61) 
(98) 


o.90ll°;°f 7 

o.908±8:8i| 


(38) 
(37) 


0.523 (180) 
0.530 (177) 


0.377 (199) 
0.382 (199) 


0.835 (3) 
0.826 (4) 



Note. — Numbers in the parentheses represent the number of occurrences that give larger t values than the 
measured value among the two hundred SILC400 simulations. 
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Fig. 16. — Quadrupole and octopole powers (<5T ; 2 ) and angu- 
lar separations between the quadrupole and the octopole (623) 
measured from the ILC maps produced from the WMAP data 
and the Galactic foreground maps with various angular resolution. 
Each WMAP or Galaxy map with FWHM = 1° has been further 
smoothed with a smoothing filter of 0?5 to 3?0 FWHM with steps 
of 0?5. Total resolution of ILC maps varies from 1?12 to 3?16 
FWHM. 



group index map composed of four hundred pixel-groups 
that contain pixels with similar foreground spectral prop- 
erties, and obtained a CMB map with foreground emis- 
sion effectively reduced (SILC400). 

Two hundred ILC simulations show that the residual 
foreground emission in the SILC400 is very small in am- 
plitude, which is at 10-30 /iK levels at high Galactic lat- 
itude regions (Fig. Et). On high latitude regions where 
the foreground emission is very weak at high frequency 
bands (V and W), the spectral index is inaccurately mea- 
sured due to the MEM reconstruction noise in the Galac- 
tic foreground maps. If ILC regions are defined based 
on the inaccurate spectral index information, the per- 
formance of foreground removal with the ILC method 



becomes low. To reduce such noises, we have smoothed 
the foreground maps with a Gaussian filter whose width 
varies depending on the signal-to-noise level of intensities 
in the MEM Galaxy maps. 

Such residual foreground features do not affect sig- 
nificantly the statistics of large-scale modes of CMB 
anisotropy. For example, when the average residual emis- 
sion map (Fig. [8J1) is added to or subtracted from the 
SILC400 map, the estimated quadrupole power and an- 
gular separation are 5T% — 288 /iK 2 and #23 = 13° for 
the addition map and 6T$ = 264 /iK 2 and 9 23 = 11° 
for the subtraction map. Both are within the 68% con- 
fidence limits of the SILC400 results. This implies that 
the residual foreground emission in our analysis map is 
not statistically important to the large-scale modes of 
CMB anisotropy. 

We have also shown that the SILC400 recovers the 
true low ^-mode powers, angular separation between the 
quadrupole and the octopole, and the planarity parame- 
ter (Figs. El-DHEIl andHU). Our SILC400 correctly re- 
constructs the spherical harmonic coefficients with min- 
imal biases. Our ILC method is insensitive to the num- 
ber of pixel-groups of common spectral index and on the 
smoothing scale as demonstrated in Figures El an( l H6I 

The quadrupole and octopole powers measured from 
SILC400 are 244l^ 3 fjK 2 and 859±fo 9 /iK 2 (68% C.L.), 
respectively. According to the SILC400 simulations, 
the minimum variance ILC method tends to underes- 
timate the low Z-mode powers. Removing the effect of 
bias due to the residual Galactic emission on the spher- 
ical harmonic coefficients of the SILC400, we obtain 



8T4 



276 +94 

z '°-126 



^K 2 and STi 



952±j$ (iK* 



We con- 



firm that the CMB quadrupole power is still lower than 
the theoretical value of the concordance ACDM model, 
with p = 5.7%. Our estimate is consistent with that of 
WILC3YR, and is located between those measured from 
the TCM3YR and LILC. 

The quadrupole power of the SILC400 is consistent 
with the previous meas urements from th e high lati- 
tude part of CMB data. lEfstathioul l)2004|) applied the 
maximum likelihood analysis method to the foreground- 
reduced CMB maps, and measured the quadrupole pow- 
ers, 223 /iK 2 for WIL C1YR and 250 ^K 2 f or TCM for 
the Kp2 sky coverage. iBielewicz et al.l l|2004D applied the 
power equalization filter to the WMAP data with Kp2 
mask, and obtained quadrupole powers 248 ± 45 /iK 2 for 
WILC1YR and 279 ± 49 fiK 2 for TCM. 
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The angular separation between the quadrupole and 
the octopole and the their planar properties have 
also been investigated. We have confirmed that the 
quadrupole and the octopole are aligned with high pla- 
narity. The probabilities of observing such anomalies 
from the bias-corrected SILC400 map are 4.3% for angu- 
lar separation #23 and over 18% for planarity parameters 
(see Tables 0] and |5j). The observed angular separation is 
marginally statistically significant. 

As observed in Figure ^1 the power spectra of the 
CMB and the Galactic emission maps show zigzag pat- 
terns at I = 2-8, with the zigzag directions opposite to 
each other, which strongly implies that the large scale 
CMB signal estimated from SILC400 and other analyses 
may be contaminated by the Galactic foreground emis- 
sion and is anti-correlated with the Galactic signal (see 
also Tables Hand©- First, we have investigated whether 
the observed zigzag pattern at I = 2-8 in the CMB power 
spectrum is statistically significant or not. Among the 
200 true input CMB signal maps, we have found 16 cases 
with the CMB zigzag pattern, 10 cases with the Galaxy 
zigzag direction, and totally 26 cases (13%). 

Secondly, we have measured cross-correlation at zero 
lag C(0) = ATs(ni)T G (ni)/N p between the SILC400 
(bias-corrected) containing only I = 2-8 modes and the 
K-band MEM Galactic foreground map smoothed by the 
FWHM = 7° Gaussian filter. Both maps have been de- 
graded to Nside = 32 before cross-correlation measure- 
ment. The C(0) values estimated from pixels at \b\ < 20° 
are —0.008 and —0.020 mK 2 for northern and southern 
hemispheres, respectively. We see that the CMB signal 
on the southern hemisphere has strong anti-correlation 
with the Galactic emission. We have calculated the cross- 
correlation functions between the 200 true input CMB 
signal maps and the Galactic emission map. The median 
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values together with 68% confidence limits for C(0) are 
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and southern hemispheres, respectively. The estimated 
probabilities that the large scale modes of true input 
CMB signal map have anti-correlation with the Galactic 
emission map as low as the measured values by chance 
are 36.5% and 19.5% for northern and southern hemi- 
spheres, respectively. The large scale modes of SILC400 
on the southern hemisphere also shows anti-correlation 
with the Galactic emission even at \b\ > 20°. However, 
whether such anti-correlation occurs due to the residual 
Galactic emission or by chance is not clear and further 
investigation is needed. 

In this study, to measure the spectral indices of the 
foreground, we have used the MEM-derived Galactic 
foreground maps that contain some level of reconstruc- 
tion errors and fail to model the Galaxy on the Galac- 
tic plane region. For the best reconstruction of CMB 
anisotropy through the ILC method, it is essential to 
use the Galactic emission maps that model the Galaxy 
realistically for precise measurement of the foreground 
spectral indices. 
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